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ABSTRACT 

The  plane  strain  vibration  frequencies  of  an  infinitely  long  hollow  cylinder 
are  calculated  exactly  with  the  aid  of  a  high-speed  electronic  computer,  for  a 
range  of  wall  thicknesses  and  azimuthal  node  numbers,  and  for  a  variety  of 
boundary  conditions.  These  latter  cover  the  cases:  outer  radial  surface  free, 
or  supported,  or  clamped;  inner  surface  free,  or  supported,  or  matched  to  a 
gas  filling  the  inner  cavity.  In  the  last  case,  we  calculate  the  attenuation  of 
sound  in  the  cylinder  wall,  in  addition  to  the  eigenfrequencies.  Our  computa¬ 
tions  of  mode  frequencies  disclose  that:  (i)  the  shear -compressional  wave 
coupling  in  the  solid  is  very  weak  (except  for  low  frequency);  (ii)  for  the  gas- 
filled  cylinder,  the  modes  are  very  close  either  to  the  gas  quasi-modes  (gas 
in  a  rigid -walled  cylinder),  or  to  the  solid  qv*si-modc3  (shell  with  free  inner 
surface),  depending  on  the  geometry. 

I.  INTRODUCTION 

A  fundamental  problem  of  the  classical  theory  of  elasticity  is  the  determina 
tion  of  the  normal  modes  of  elastic  vibration  of  a  hollow  cylinder.  In  principle, 
this  problem  is  solved  by  writing  down  the  solution  of  the  elastic  wave  equation 
and  then  determining  the  mode  frequencies  by  the  requirement  that  an  appro¬ 
priate  set  of  boundary  conditions  be  fulfilled.  In  practice,  the  problem  is 
unfortunately  not  so  simple. 

The  major  difficulty  is  that  the  determination  of  the  characteristic  frequen¬ 
cies  generally  requires  an  enormous  amount  of  numerical  work,  particularly 
because  the  secular  equation  involves  cylinder  functions  whose  arguments 
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depend  on  the  frequency  to  be  computed.  Since  the  cylinder  problem  possesses 
a  multitude  of  independent  parameters  (cylinder  wall  thickness,  elastic  con¬ 
stants,  and  wave  numbers)  and  since  a  considerable  number  of  sets  of  boundary 
conditions  must  be  investigated,  it  is  clear  that  a  complete  solution  of  the 
problem  would  involve  a  formidable  program  of  numerical  computation. 

Of  course,  in  simple  special  cases,  it  is  possible  to  obtain  the  modes  for 
any  wall  thickness  from  the  exact  theory,  without  undue  labor.  Such  calcula¬ 
tions  previously  have  been  done  for  the  purely  radial  dilational  modes  of  free, 

infinitely  long  cyh  ders,  *  and  for  the  purely  axial  shear  modes  of  free  and 

2 

clamped  infinite  cylinders. 

Now  it  is  feasible,  with  electronic  computers,  to  calculate  exactly  the 

vibration  modes  for  thick-walled  cylinders  ir.  other  cases  as  well.  Already 

3 

Herrmann  and  Mirsky  have  presented  such  calculations  for  axially  symmetric 

longitudinal  waves  in  free,  infinitely  long  cylinders.  They  gave  numerical 

4 

results  for  the  lowest  mode  and  a  range  of  wall  thicknesses.  Greenspon 

S.  Field,  Can.  J.  Research  A 17,  141  (1939);  J.  A.  McFadden, 

J.  Ac.oust.  Soc.  Am.  26,  714  (1954);  A.  E.  Bassett,  Proc.  London  Math. 

Soc.  21,  53  (1889). 

2 

J.  H.  Baltrukonis  and  W.  G.  Gottenberg,  J.  Acoust.  Soc.  Am.  31^, 

734  (1959). 
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G.  Herrmann  and  I.  Mirsky,  J.  Appl.  Mech. ,  Trans.  Am.  Soc. 

Mech.  Engrs.  78,  563  (1956),  79,  1  (1957). 
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J.  E.  Greenspon,  J.  Acoust.  Soc.  Am.  31,  1682  (1959),  J.  Aero.  Space 


Sci,  27,  37  (1960). 
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calculated  the  lowest  few  flexural  modes  for  free,  infinitely  long  cylinders. 

In  addition,  he  has  treated  the  flexural  modes  of  the  finite  length  cylinder, 

again  with  free  sides,  but  with  ends  either  free  or  freely-supported.  More 

recently,  Gazis  has  published  extensive  calculations  for  the  infinite  hollow 

cylinder  which  is  stress-free  at  both  inner  and  outer  radial  surfaces  -  first, 

for  plane-strain  vibrations,  later  for  the  general  case  involving  axial  motion 
5 

also.  He  considered  a  range  of  geometries  and  axial  wavelengths  and  several 
values  of  the  azimuthal  quantum  number  for  the  lowest  dozen  modes  or  so. 

The  purpose  of  the  present  paper  is  to  calculate  exactly  the  eigenfrequencies 
of  the  infinite  length  thick-walled  cylinder  for  boundary  conditions  not  treated 
in  the  above  work  (refs.  3,  4,  5),  i.e. ,  for  other  than  stress-free  radial  sur¬ 
faces.  In  these  computations,  we  shall  consider  plane -strain  vibrations  in  the 
infinite  hollow  cylinder  for  a  range  of  wall  thicknesses.  For  these  modes, 
however,  we  shall  extend  Gazis'  work  as  follows. 

First  of  all,  the  infinite  hollow  cylinder  will  be  treated  for  a  number  of  new 
combinations  of  boundary  conditions  of  the  usual  type.  Thus,  we  shall  consider 
the  outer  surface  to  be  either  clamped,  r  supported,  or  free,  while  the  inner 

surface  is  either  supported,  or  free.  These  investigations  lead  to  the  discovery 

6 

that  the  thickness  eigenmodes  are  generally  either  almost  purely  radial  dila- 
tational  or  almost  purely  azimuthal  shear.  The  frequencies  of  these  modes  are 

^D.  C.  Gazis,  J.  Acovst.  Soc.  Am.  JO,  786  (1958),  _31,  568,  573  (1959). 
®These  are  modes  with  at  least  one  displacement  node  or  antinode 
occurring  across  the  thickness  of  the  cylinder  wall. 
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in  fact  very  nearly  the  same  as  for  the  corresponding  pure  modes.  This  indi¬ 
cates  that  the  shear-compression  coupling  is  quite  weak  for  such  modes. 

Secondly,  we  shall  considei  tne  very  interesting  situation  of  a  gas  filling 

7  8 

the  inner  part  of  the  cylinder,  1  while  the  outer  surface  is  again  either 
clamped,  or  supported,  or  free.  In  this  case  we  are  dealing  with  the  acoustic 
modes  of  a  two-phase  system,  where  no  matter  how  tenuous  the  gas  inside  the 
cylinder  may  be,  the  modes  of  the  system  are  often  drastically  different  from 
those  calculated  by  neglecting  the  presence  of  the  gas.  We  find  that,  depending 
on  the  ratio  of  inner  to  outer  radius,  the  eigenmodes  are  very  closely  similar 
either  to  the  modes  of  a  shell  with  free  inner  surface  (solid  quasi-mode)  or  to 
the  modes  of  a  gas  in  a  rigid -walled  cylinder  (gas  quasi-mode). 

7 

Some  special  examples  of  this  case  have  been  given  previously  in  a 
theoretical  treatment  of  the  problem  of  combustion  instability  in  a  rocket  motor 
charged  with  an  inside -burning  tubular  grain.  See,  (a)  F.  T.  McClure,  R.  W. 
Hart,  and  J.  F.  Bird,  J.  Appl.  Phys.  (scheduled  for  publication  in  May  1960), 
and  (b)  F.  T.  McClure,  R.  W.  Hart,  and  J.  F.  Bird,  Proceedings  of  the 
Princeton-A. R. S.  Solid  Propellant  Research  Symposium  (Jan.  28,  29,  1960). 

8 

Other  such  two-phase  systems  have  been  treated  previously  by  C. 
S&lceanu  and  M.  Z&g&nescu,  Compt.  rend.  247,  812  (1958)  and  F.  I.  N. 
Niordson,  Trans.  Roy.  List.  Tech.  Stockholm,  no.  73  (1953)  for  very  simple 
cases,  and  by  R.  D.  Fay,  J.  Acoust.  Soc.  Am.  24,  459  (1952)  with  a  vitiating 
error.  Also  see  M.  A,  Biot,  J.  Appl.  Phys.  23,  997  (1952). 
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Finally,  as  a  rather  novel  application  of  elasticity  theory,  we  compute 
exactly  the  attenuation  of  sound  in  the  solid  wall  for  the  gas-filled  cylinder, 
by  only  slightly  extending  the  eigenfrequency  computation  program. 


H.  DISPLACE  ME  NT  AND  STRESS  IN  CYLINDER 

Here,  we  shall  merely  display  the  form  of  the  solution,  since  the  methods 

9 

of  derivation  are  standard  procedure  and  readily  available  elsewhere.  In 
the  notation  of  ref.  7,  we  have  for  the  displacement  vector  J3,,  and  stress 
tensor,  P  ,  the  following  components: 

.  iwt  ,  .  r.  Tt  A  /I-  _\ 


Sr  =  +  ie  cos(m<p)  j  A^^r)  +  A2Ym(kcr) 


+  ITT  P¥m<V>  +  V»<V>J 

s 


•  pia,t  r 

-  -P1--  sin(m</7)  VA,J  (krl+A.Y  (kr) 

kcr  '  '  )  1  m'  c  •  2  mv  c  ’ 


[B,  J*  (k  r)  +  B0Y’  (k  r)l 
L  1  m'  s  2  m  s  ,J 


9  iwt  r 

Prr  = - E -  cos(m<?)  J  A1jr(kcr)  +  A2yr(kcr) 


Bl3l(k3r)  +  B2yi(k8r) 


%  =  +  “7 -  sin<ra  &  *  A2yl(V> 

+  Bl^  <ksr)  +  B2y«,  <ks">J^  <4> 


P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics  (McGraw- 
Hill  Book  Company.  Inc.,  New  York,  1953),  pp.  142,  1840. 
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where 


jj(x)  =  m 


Jm(x) 

m  Tt  /  » 

—  -  Vx) 


jrw  *  JaW  +  x  Jm(x) 


i  + 


sr  <V>2  -  “2 
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and  the  y  functions  similar  but  with  J - VY  .  Also, we  have  used, 

m  '  m 

2  «  2 


u2  _  _ 

*©  t  '  V;t" 
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.2  &> 
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i 
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(8a) 


with  the  familiar  relationships 


where  X  and  p  are  the  Lame  elastic  moduli,  p  is  the  density  of  the  solid 
cylinder,  and  m  is  the  (integer)  angular  wave  number.  The  determination  of 
the  amplitudes  A^,  Ag,  3^,  Bg  by  applying  the  boundary  conditions  accom¬ 
plishes  the  solution  to  the  problem. 


m.  THE  BOUNDARY  CONDITIONS 

In  this  section  we  specify  the  various  sets  of  boundary  conditions  that  we 
shall  apply  to  the  displacement  and  stress  functions,  Eqs.  (1-4)  9  the  inner 
and  outer  boundaries  of  the  cylinder,  i.  e. ,  at  r  -  a,  r  =  b.  (Since  the  cylinder 
length  is  infinite,  we  apply  no  boundary  conditions  on  planes  perpendicular  to 
the  x-axU.) 
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1.  Conditions  on  Outer  Cylindrical  Surface 
On  this  surface  r  =  b,  we  shall  apply  any  one  of  the  three  sets  of  condi - 
§  tions 


P  =P„  =0, 

at  r  =  b, 

(9a) 

rr  rxp  ’ 

S  =  P  =0, 
r  T(p 

at  r  =  b, 

(9b) 

ii 

w 

it 

© 

at  r  =  b, 

(9c) 

corresponding  to  free,  supported,  or  clamped  ouside  surface,  respectively. 


2.  Conditions  on  Inner  Cylindrical  Surface 
The  inner  surface  (r  =  a)  of  the  hollow  cylinder  we  shall  consider  to  be 
either  free 


Prr  =  %  =  0,  at  r  =  a  ,  (10a) 

or  supported 

Sr  =  %  =  °>  at  r  =  a>  <10b) 

or  to  match  in  both  stress  and  displacement  to  a  gas  filling  the  inside  of  the 
cylinder.  Assuming  the  ga3  exerts  no  tractive  force  on  the  solid  surface,  this 
last  case  requires 


rr 


+  ia)ZS  =P  =  0,at  r  =  a 
g  r  rep  » 


T(p 


(10c) 


where  Zg  is  the  acoustic  impedance  for  the  gas  at  the  boundary  wall, 

Z  2.  [p/u]  _  ,  with  p,  u  being  acoustic  pressure  and  radial  particle  velocity 

g  r-a 

in  the  gas.  With  acoustic  gas  pressure  of  the  form  p^1  cos  (m^JJ^k^r),  we 


have  from  u  = 


— - —  *  that  the  impedance  is 

Pgo>  3r 
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where  p  is  gas  density,  and  k  =  w/c  ,  with  c  being  sound  velocity  in  gas. 

8  STS 

We  note  that  Eqs.  (10a),  (10b)  are  special  cases  of  (10c),  for  Z  ■ -*0  and 

o 

Z  ~>co,  respectively. 

B 

When  we  apply  these  various  sets  of  boundary  conditions  to  the  elastic 
wave  described  by  the  displacement  and  stress  of  Eqs.  (1)  »  (4)  ,  we  will 
find  in  each  case  that  four  simultaneous  equations  which  are  linear  and 
homogeneous  in  the  constants  A^,  Ag,  Bj,  B2,  need  to  be  satisfied.  For  a 
non-trivial  solution,  the  determinant  formed  from  the  coefficients  of  the  A's 
and  B's  must  be  zero.  This  condition  determines  the  eigenfrequencies  of  the 
system,  so  that  in  the  following  we  are  to  be  concerned  with  finding  the  fre¬ 
quency  roots  of  this  secular  determinant  for  each  set  of  boundary  conditions. 

IV.  DIMENSIONS  AND  PHYSICAL  CONSTANTS  OF  THE  SYSTEM 

The  secular  equation  for  the  oscillation  frequencies  (f  =  o>/2ir)  will  in 

general  involve  the  dimensions  a,  b,  of  the  cylinder,  and  the  density  and  elastic 

constants  of  both  the  cylinder  and  the  gas  inside  the  cylinder.  For  the  empty 

cylinder,  these  quantities  fall  into  three  dimensionless  groups,  a/b,  wb/c  ,  and 

s 

a  (Poisson’s  ratio);  for  the  gas  filled  cylinder,  two  additional  parameter  groups 
appear,  namely  Pg/p,  and  c g/cs*  Clearly,  a  complete  investigation  over  the 
whole  parameter  field  would  require  an  undesirably  large  effort.  In  this  modest 
investigation,  we  have  confined  our  attention  to  a  study  of  the  dependence  on 
only  a/b  and  wb/cg,  while  selecting  *he  constants  o  =  0. 38,  p ^/p  -  1/160,  and 
Cg/cg  =  10/3.08  as  representative  values  of  the  remaining  quantities.  Rather 
than  present  the  results  in  terms  of  the  two  dimensionless  parameters  which 
we  have  considered,  however,  we  elect  to  display  the  mode  frequencies  as  a 


s 


function  of  inner  radius  for  a  selected  value  of  b,  and  c  ,  because  we  believe 

that  this  presentation  facilitates  a  more  concrete  feeling  for  the  nature  of  the 

4 

results.  We  shall  use  b  =  5. 75  cm,  cg  =  3.08  x  10  cm/ sec.  (These  values 

a  2  9  2k 

correspond  to  A  =  4. 313  x  10°  dyne/cm  and  =  1.  52  x  10  dyne/cm 

Finally,  in  our  illustrative  calculations  of  sound  attenuation  in  the  solid  wall, 

we  shall  use  the  dilatational  and  shear  viscosities  A*  =  2C0  poise,  77  =  200  poise, 

7 

respectively. 

V.  EXACT  NUMERICAL  EIGENFREQUENCIES 
1.  Clamped,  Supported  and  Free  Boundaries 
Here  we  present  our  results  for  the  six  sets  of  boundary  conditions  speci¬ 
fied  by  Eqs,.  (9),  (Ida),  (10b),  viz. ,  outer  surface  either  clamped  or  supported 
or  free,  inner  surface  either  supported  or  free. 

The  characteristic  frequencies  were  calculated  on  the  ’’Univac -Scientific” 
for  several  values  of  m .  Roots  were  sought  for  frequencies  up  to  25,000  cps , 
corresponding  to  wb/c  -C  29. 4.  The  method  of  computation  was  to  select  a 
value  of  a  and  then  calculate  the  deter minantal  roots  by  an  interval-halving 
technique.  Doing  this  for  about  10  values  of  a  from  2. 87  to  5. 75  cm  was  suffi¬ 
cient  to  give  the  mode-maps  for  these  cases,  since  the  curves  are  smooth. 

These  results  are  shown  in  Figs.  1  to  6,  as  graphs  of  frequency  vs  inner 
radius/ outer  radius. 

Fig.  I  shows  the  plane -strain  modes  of  the  completely  free  cylinder  (cf. 

refs.  4, 5).  Those  modes  such  that  the  frequency  f  becomes  infinite  as 

4 

a-*b  are  the  thickness  modes  defined  in  Sec.  I.  We  follow  Gazis  in  calling 
the  remaining  pair  of  modes  ’’ring  modes'*.  Figures  2-6  illustrate  how  different 
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Fig.  5  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  m  -  0,  1  AND  2 
MODES  OF  A  SOLID  CYLINDRICAL  SHELL  WITH  OUTER  SURFACE 
CLAMPED  AND  INNER  SURFACE  FREE 
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sets  of  boundary  conditions  affect  both  thickness  and  ring  modes. 

Consider  first  the  ring  modes.  The  lower  of  these  in  Fig.  1  is  a  flexural 

mode,  which  has  non-  zero  frequency  only  for  m^2,  while  the  upper  is  an 

4 

extensional  (,4breathing")  mode.  For  each  of  Figs.  2-6,  the  extensional 
freedom  is  removed  (8^  =  0  at  r  =  a  and/or  r  =  b)  and  the  ring-extensional  mode 
disappears.  In  Figs.  5  and  6  the  ring-flexural  mode  is  also  absent,  since  the 
flexural  freedom  has  been  removed  (S^  =  0  at  r  =  b).  Note  that  for  Figs.  2,  3, 
and  4  the  ring-flexural  mode  has  zero  frequency  only  for  m  =  0. 

The  thickness  modes  also  exhibit  a  regular  behavior  as  boundary  conditions 
are  changed.  One  group  of  these  modes  (labeled  "S"  in  Figs.  1-6)  has  pre- 
dominantly  shear  wave  characteristics,  while  the  remaining  thickness  modes 
(labeled  "U'  in  Figs.  1-6)  form  a  group  with  predominantly  dilatational  wave 
character.  This  is  particularly  evident  from  the  manner  in  which  the  frequency 
curves  for  each  group  are  affected  by  the  boundary  conditions.  Thus  the  S-mode 
curves  for  free  and  for  supported  boundaries  are  very  nearly  identical,  while 
both  of  these  are  vastly  different  from  the  S  curves  for  clamped  boundaries 
(see  the  two  sets  of  figures,  1-4  and  5-6).  On  the  other  hand,  the  D-mode 
curves  for  supported  and  for  clamped  boundaries  are  strikingly  alike  while 
both  of  these  differ  widely  from  the  D  curves  for  free  boundaries***  (see  the  two 

***The  similarity  of  the  D  curves  of  Fig.  1  to  those  of  Figs.  3,6, 
despite  the  difference  in  boundary  conditions,  is  due  to  the  fact  that  in  each  of 

these  cases  the  boundary  conditions  imply  n  =  b  -  a,  n  =  1,  2,  3 _ ,  where 

X  is  the  (ill-defined)  radial  wavelength. 
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sets  of  figures,  3,  6  and  2, 4, 5). 

These  conclusions  are  fortified  by  an  inspection  of  the  machine  calcula¬ 
tions  for  each  mode.  For  this  discloses  that  indeed  the  S-modes  have  displace¬ 
ment  and  stress  fields  largely  shear  in  character  (terms  involving  Bp  B ^ 
dominate  in  Eqs.  (l)-(4^  while  the  D-modes  have  displacement,  and  stress 
largely  dilatational  in  character  (terms  with  A^Ag  dominate  in  Eqs.  (1)  -  (4jjl 

It  appears  then  that  the  eigenfrequencies  of  the  S-modes  and  D-modes  are 
astonishingly  close  to  those  of  pure  shear  and  pure  dilatational  waves,  respec¬ 
tively  (except  for  the  lowest  modes  of  very  thick  shells).  To  a  lesser  degree, 
the  stress  and  displacement  fields  are  similar  to  those  of  pure  modes.  Thus, 
it  is  clear  that  the  shear  and  dilatational  waves  are  only  very  weakly  coupled 
by  our  boundary  conditions,  and  that  the  "frequency  mixing”  is  extremely  small. 
These  statements  are,  of  course,  exactly  true  for  axially  symmetric  modes 
(m  =  0  in  Eqs.\(l)  r  (.4),  also  cf*  ref.  4). 
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2.  Gas-Filled  Cylinder 

We  now  turn  to  the  mode  maps  for  the  three  sets  of  boundary  conditions 
given  by  Eqs.  (9)  with  (10c),  corresponding  to  a  cylinder  filled  with  a  '’stiff* 
gas,  while  its  outer  surface  is  either  free,  or  supported,  or  clamped. 

The  numerical  work  was  considerably  more  complicated  here  than  for  the 
cases  described  previously  (Sec.  V,  1).  Thic  is  because  the  fvsa  curves  for 
the  gas-filled  cylinder  are  somewhat  jagged,  aud  at  some  points  in  the  f -a 
plane  two  modes  come  close  together.  The  method  of  computation  adopted  was 
to  trace  out  each  mode  curve  separately  as  a  function  of  ja,  by  a  slope -projection 
technique.  Where  two  modes  approach  each  other  closely,  accidental  jumping 
from  one  to  the  other  was  avoided  by  collapsing  the  interval  until  Z  (k^a) 
varied  smoothly.  We  have  done  these  calculations  for  several  values  of  m  . 

In  ref.  7,  the  mode  maps  for  m  =  0,  clamped  outer  surface,  and  m=  2, 
clamped, supported  or  free  outer  surface,  have  already  been  presented.  Here 
we  display  in  Figs.  7-12  the  mode  maps  for  m  =  0, 1, 3  not  heretofore 
published. 

The  axially  symmetric  (m  =  0),  or  pure  radial  modes  are  shown  in  Figs. 

7,8  for  free  and  supported  outer  surface,  respectively.  We  have  plotted  only 
the  dilatational  modes,  since  in  this  case,  the  shear  modes  are  not  coupled  to 
the  gas  and  thus  are  the  same  as  in  Figs.  1-6.  As  might  be  expected,  the 
mode  frequencies  correspond  rather  closely  to  either  the  solid  or  gas  quasi - 
modes,  i.e. ,  the  modes  of  a  solid  shell  with  free  inner  surface,  or  the  modes 
of  a  gas  in  a  rigid  walled  cylinder,  respectively.  These  quasi -modes,  where 
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Fig.  7  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  RADIAL  MODES  (m=0'‘ 
OF  THE  GAS-FILLED  CYLINDER  WITH  FREE  OUTER  RADIAL  SURFACE 
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Fig.  8  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  RADIAL  MODES 
(m  -  0)  OF  THE  GAS-FILLED  CYLINDER  WITH  SUPPORTED 
OUTER  RADIAL  SURFACE 
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Fig.  9  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  FIRST  AZIMUTHAL 
MODES  (rc  -  1)  OF  THE  GAS-FILLED  CYLINDER  WITH  FREE 
OUTER  RADIAL  SURFACE 
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Fig.  10  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  FIRST  AZIMUTHAL 
MODES  (m  -  1)  OF  THE  GAS-FILLED  CYLINDER  WITH  SUPPORTED 

OUTER  RADIAL  SURFACE 
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V ig.  iJL  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  FIRST  AZIMUTHAL 
JfcODES  (m  «  1)  OF  THE  GAS-FILLED  CYLINDER  WITH  CLAMPED 

OUTER  RADIAL  SURFACE 
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different  from  the  actual  modes,  are  indicated,  for  example,  by  the  dashed 
lines  in  Fig.  8.  Since  the  interface  boundary  conditions  for  the  quasi-modes 
are  mutually  exclusive,  the  modes  of  the  two-phase  system  avoid  the  inter¬ 
sections  of  dashed  lines,  and  alter  their  character  in  such  neighborhoods  from 
quasi-modes  of  one  medium  to  those  of  the  other. 

Figs.  9,  10,  11  show  the  first  azimuthal  modes  {m  =  1)  for  free,  supported 
and  clamped  outer  surface,  respectively.  Here  the  modes  of  the  two-phase 
system  are  not  individually  pure  shear  or  pure  dilatationai,  as  they  are  when 
m  =  0,  but  combinations  of  both,  so  that  the  mode -maps  appear  rather  compli¬ 
cated.  We  note  that  the  '’splitting*'  of  modes  at  gas  quasi-mode  and  solid  quasi - 
mode  intersections  is  influenced  not  only  by  the  mismatch  in  acoustic  proper¬ 
ties  of  the  solid  and  the  gas,  but  also  by  the  relative  amounts  of  shear  and 
dilatationai  motion  in  the  solid.  For  example,  compare  Fig.  11  with  Fig.  5, 
where  we  see  that  where  a  gas  quasi-mode  intersects  a  type  "D"  solid  quasi¬ 
mode,  the  splitting  is  much  larger  than  where  it  intersects  a  type  "S"  solid 
quasi-mode.  This  may  be  presumed  to  occur  because  the  gas  exerts  no  shear 
on  the  solid,  and  couples  directly  only  to  the  dilatationai  motion  -  and,  as  we 

saw  in  Sec.  V.  1,  the  "S-modes”  have  little  dilatationai  character,  while  the 

11 

"D-modes"  are  largely  dilatationai. 

Finally,  Fig.  12  shows  the  third  azimuthal  (m  =  3)  modes  of  the  cylinder 
with  supported  outer  radial  surface. 

**For  further  qualitative  discussion  of  the  modes  of  the  two-phase  system, 


see  ref.  7b.*. . 
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VI.  ATTENUATION  OF  SOUND  IN.  CYLINDER  WALL 


We  shall  now  calculate  the  rate  of  dissipation  of  elastic  wave  energy  in 

the  wall  of  the  gas-filled  cylinder  (cf.  ref.  7a).  In  a  stationary  state,  this 

is  equal  to  the  rate  at  which  energy  flows  into  the  cylinder  wall  across  the 

12 

inner  surface  (r  =  a),  so  that,  in  the  acoustic  approximation,  the  dissipation 
rate  per  unit  length  of  solid  wall  is 


(kga)  |  Po 


Real  part  of 


r  =  a 


where  pQ  is  the  gas  pressure  amplitude  (Sec.  ID.  2).  The  real  part  of  the 
boundary  admittance  is  given  to  first  order  in  the  coefficients  of  expansive 
friction  (X)  and  shear  viscosity  (77)  by^a 


We  note  that,  for  the  physical  constants  used  in  this  work  (Sec.  IV),  the  first 
order  expression  (12)  should  be  accurate  to  a  few  percent  up  to  frequencies 
/v  100  kc. 

The  derivatives  in  Eqs.  (13)  were  calculated  during  the  process  of  finding 
the  characteristic  frequencies  of  the  gas-filled  cylinder  (Sec.  V.  2).  Once  the 
eigenfrequency  for  a  given  a  was  determined,  A  (or  /x)  was  varied  slightly 
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The  energy  flux  at  r  =  b  is  zero  for  all  boundary  conditions  treated  in 
this  paper  (see  Sec.  III). 


*x* 
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and,  for  the  same  a  and  the  same  frequency,  the  coefficients  Ag,  Bj, 

Bg  were  recalculated  as  minors  of  the  secular  determinant.  The  deriva¬ 
tives  are  then  obtained  from  Eqs.  (1),  (3),  and  the  real  part  of  the  boundary 
admittance  from  Eq.  (12).  The  dilatational  part  of  the  damping  (A  -  200  poise, 

7]  -  0)  and  the  total  damping  (A*  =  r\  -  200  poise)  are  shown  separately  in  the 
figures. 

Some  attenuation  calculations  for  the  m  =  2  modes  have  already  been  pre¬ 
sented  in  ref.  7a.  Here  we  display  the  results  for  the  m  =  0  and  m  =  1  modes 
of  the  cylinder  with  supported  outer  surface  (Figs.  13  and  14).  The  admit¬ 
tances  are  plotted  as  functions  of  a/b  along  those  portions  of  the  modes  which 
are  closest  to  the  lowest  gas  quasi-mode  (cf.  Figs.  8  and  10),  since  they  are 
very  large  except  very  close  to  the  gas  quasi  *  mode.  In  Fig.  15,  we  have  plotted 
the  admittance  for  a  first  azimuthal  ring -mode  (cf.  lowest  curve  in  Fig.  10)  for 
the  supported  boundary  at  r  =  b.  Note  that  the  damping  is  about  10,000  times 
larger  for  the  ring-mode  than  for  the  thickness  modes. 
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Fig.  13  ATTENUATION  OF  SOUND  IN  THE  SOLID  WALL  OF  THE  GAS 
FILLED  CYLINDER  WITH  SUPPORTED  OUTER  SURFACE 
FOR  RADIAL  MODES  (m  -  0)  LYING  NEAR  TO  THE 
LOWEST  GAS  QUASI- MODE 


F  SPECIFIC  ACOUSTIC  ADMITTANCE  (royl) 


OL — 
0500 


l  *  \ 

t  )  ' 


0.625 


0.750 


0.375 


I.OCO 


RATIO  OF  INNER  RADIUS  TO  OUTER‘-r- 

o 


Fig.  14  ATTENUATION  OF  SOUND  IN  THE  SOLID  WALL  OF  TEE  GAS- 
FILLED  CYLINDER  WITH  SUPPORTED  OUTER  SURFACE 
FOR  FIRST  AZIMUTHAL  MODES  (m  =  1)  LYING 
NEAR  TO  THE  LOWEST  GAS  QUASI-MODE 
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VII.  CONCLUSION 

In  this  paper  we  have  been  concerned  primarily  with  exact  calculation  of 
the  plane-strain  eigenfrequencies  of  an  infinitely  long,  hollow  cylinder  of 
arbitrary  wall  thickness.  We  have  considered  free,  supported  and  clamped 
radial  surfaces  in  heretofore  unstudied  combinations.  In  addition,  we  have 
presented  new  mode-maps  (frequency  vs  geometry)  for  a  gas-filled  cylinder. 
Finally,  we  have  evaluated  the  attenuation,  of  sound  energy  in  the  thick 
cylindrical  shell. 

From  a  study  of  our  results  for  the  diverse  boundary  conditions  it  has 
appeared  that  the  thickness  eigenmodes  may  be  regarded  as  consisting  of  two 
distinct  classes,  the  first  being  almost  purely  radial  dilatational  waves,  (D 
class),  the  second  almost  purely  azimuthal  shear  waves,  (S  class).  The 
sharpness  of  this  dichotomy  evidences  very  weak  shear ~compressional 
coupling,  so  that  the  thickness  modes  should  be  easily  obtainable  from  the  pure 
modes  by  perturbation  theory.  The  development  of  such  an  approximate  theory 
will  be  the  topic  of  a  subsequent  paper  by  one  of  the  authors. 

The  calculations  for  the  gas -filled  cylinder  allowed  that  the  mode  frequen¬ 
cies  are  very  close  either  to  those  of  the  gas  quasi-modes  or  to  those  of  the 
solid  quasi-modes,  depending  on  the  radial  dimensions  of  the  cylinder  and  the 
number  of  nodes  in  the  radial  direction. 
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